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ABSTRACT 

We numerically examine centrifugally supported shock waves in 2D rotating accretion 
flows around a stellar-mass (WMq) and a supermassive (1O 6 M0) black holes over a 
wide range of input accretion rates of 10 7 M/Mg ^ 1CP 4 . The resultant 2D-shocks 
are unstable with time and the luminosities show quasi-periodic oscillations (QPOs) 
with modulations of a factor of 2 - 3 and with periods of a tenth seconds to several 
hours, depending on the black hole masses. The shock oscillation model may explain 
the intermediate frequency QPOs with 1-10 Hz observed in the stellar-mass black 
hole candidates and also suggest the existence of QPOs with the period of hours in 
AGNs. When the accretion rate M is low, the luminosity increases in proportion to the 
accretion rate. However, when M exceeds greatly the Eddington critical rate Me, the 
luminosity is insensitive to the accretion rate and is kept constantly around ~ 3Le- 
On the other hand, the mass-outflow rate Mi oss increases in proportion to M and it 
amounts to about a few percent of the input mass-flow rate. 
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1 INTRODUCTION 

Rotating inviscid and adiabatic accretion flow around a 
black hole can have two saddle type sonic points. After the 
flow with adquate injection parameters passes through the 
outer sonic point, the centrifugal force can virtually stop the 
rotating supersonic flow, forming a standing shock close to 
the black hole. Then, the flow passes through the inner sonic 
point and again falls into the black hole sup ersonically. In 
suc h transonic problems of accretion and wind. lFukuel (| 19871 ) 
and IChakrabartil (|1989l ) firstly showed satisfactory analyt- 
ical or numerical global solutions under a full relativistic 
treatment and under a pseudo-Newtonian potential, respec- 
tively. It has been shown that these generalized accretion 
flows could be responsible for the hard and soft state tran- 
sitions or the quasi-periodic oscillations (QPOs) of the hard 
X-rays from the black hole candidates (jMolteni et al.|[l996l ; 
iRvu et al.lll997l ;l Lanzafame et al.1ll998l ). Further analyses of 
the transonic problems under modified pseudo-Newtonian 
potentials showed that the standing shocks are essential 
ingredients in multi-transonic b lack hole accretion discs 
l|Das et al.ll200ll ; [Da3l2002l . I2003T ) and that the generalized 
multi-transonic accretion model may show QPO frequencies 
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of Galacti c black hole cand idates in terms of dynamical flow 
variables <|Das et al.ll2003al lbl). 

Re cently, in rela t ive to the s hock oscillation models o f 
QPOs, lOkuda et ail (|2004l ) and IChakrabarti et al.1 (|2004l ) 
showed several numerical simulations of 2D accretion flows 
around black holes, using the Eulerian method under the 
treatment of radiation transport and the SPH method in 
the presence of cooling effects, respectively. From the power 
spectra of the luminosities, they showed that QPOs are 
found at a few Hz to hundreds of Hz for stellar-mass black 
holes with mass of WMq and oscillation time-scales of hours 
to weeks for supermassive black holes with mass of 1O 8 M0. 
In the present paper, we examine further shock oscillation 
models over a wide range of accretion rates, focusing on the 
general properties of the shock waves and the QPO frequen- 
cies in the 2D nonadiabatic accretion flows around the black 
holes. 



2 MODEL EQUATIONS 

A set of relevant equations consists of six partial differential 
equations for density, momentum, and thermal and radia- 
tion energy. These equations include the heating and cool- 
ing of gas and radiation transport. The radiation transport 
is treated in the gray, flux-limited diffusion approximation 
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l|Levermore fc Pomraninglll98ll ). We use spherical polar co- 
ordinates (r,£,<p), where r is the radial distance, £ is the 
polar angle measured from the equatorial plane of the disc, 
and ip is the azimuthal angle. The gas flow is assumed to 
be axisymmetric with respect to Z-axis (8/ dip = 0) and 
the equatorial plane . In this coordinate system, we have 
the b asic equations in the following conservative form (Kiev 
1 19891 ): 
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where p is the density, v = (v, w,v v ) are the three veloc- 
ity components, G is the gravitational constant, M* is the 
central mass, p is the gas pressure, e is the specific internal 
energy of the gas, Eo is the radiation energy density per 
unit volume, Po is the radiative stress tensor, and c is the 
speed of light. The subscript "0" denotes the value in the 
comoving frame a nd that the equa tions are correct to the 
first order of v/c jKato et al.lll998 | ). We adopt the pseudo- 
Newtonian potential (jPaczynskv fc Wiitalll980r ) in equation 
(2), where r g is the Schwarzschild radius. The force density 
/h ~ (.A - ! /c) exerted by the radiation field is given by 



ft 



: P F , 



(7) 



where k and a denote the absorption and scattering coeffi- 
cients and Fo is the radiative flux in the comoving frame. 
The quantity A describes the cooling and heating of the gas, 



A = pcn(S* - -Bo), 



(8) 



where S* is the source function. For this source function, 
we assume local thermal equilibrium S* = aT 4 , where T is 
the gas temperature and a is the radiation constant. For the 
equation of state, the gas pressure is given by the ideal gas 
law, p = RgpT/ p, where p is the mean molecular weight and 
Rg is the gas constant. The temperature T is proportional 
to the specific internal energy, e, by the relation p — (7 — 
l)pe = RopT I p, where 7 is the specific heat ratio. To close 
the system of equations, we use the flux-limited diffusion 
approximation for the radiative flux: 



and 



Ac 



P(k + a) 



grad E , 



P = E - T E 



(9) 



(10) 



where A and TEdd are the flux-limiter and the Eddington 
Tensor, respect ively, for wh ich we use the approximate for- 
mulas given in The formula fulfill the correct 



limiting conditions A — > 1/3 in the optically thick diffusion 
limit, and | F |— * c_Eo as A — > for the optically thin stream- 
ing limit. 



3 NUMERICAL METHODS 

The set of partial differential equations (l)-(6) is numeri- 
cally solved by a finite-difference method under adequate ini- 
tial and boundary conditions. The numerical schemes used 
are basically the same a s that described previously (|Klevl 
1 19891 : lOkuda et al. I ll997l ). The methods are based on an 
explicit-implicit finite difference scheme. 



3.1 Model Parameters 

We consider a stellar-mass black hole with mass M» = lOAf© 
and a supermassive black hole with M* = 1O 6 M0. To ex- 
amine the shock model, we determine the injection param- 
eters, such as the specific angular momentum, A ou t, the ra- 
dial velocity, w ut, and the sound velocity, a ou t, at an outer 
boundary radius, -Rout, whose parameters can lead to a 
shock wave close to the black hole. We search analytically 
the injection parameters through an examination of the pa- 
ramete r space (-E ou t, A ou t ) , where -E out is the t otal s pecific 
energy l|Chakrabartil 1 19891 : iMolteni et all Il994 I1999T ). The 
typical model parameters used are listed in Table 1. Here, 
the velocities and the distances are given in units of c and 
r g , respectively, and m is the input accretion rate normal- 
ized to the Eddington critical accretion rate Me(= Le/c 2 ), 
where Le is the Eddington luminosity given by 



AtyGAUc 



(11) 



and K e is the electron scattering opacity. Le and Me are 
1.5 xlO 39 erg and 1.7 x 10 18 g s -1 , respectively, for the 
stellar-mass black hole with M, — IOA^q . (f> is the subtended 
angle of the central black hole to the disc height h at r — 
Rout, that is, tan0 = (h/r) ut- The inner-boundary radius 
Ri n of the computational domain is taken to be 2r g . 



3.2 Initial and Boundary Conditions 

At the outer-disc boundary we assume a continuous inflow 
of matter with a constant accretion rate M and the injection 
flow parameters in Table 1. For the accretion rates consid- 
ered here, we have an optically thick and radiation-pressure 
dominant accretion flow at the outer boundary -Rout = 30r g . 
Then, the input gas temperature T ou t is given by 



a 4 P 4 aT\ 

^ flout J ^ o \ / out ^ r»\ Jout, 
dp 9 p 



(12) 



where P is the total pressure. Therefore, for a given a ou t, we 
have a smaller T ou t for a smaller ambient density p ou t. The 
initial conditions of the flow except the outer-disc boundary 
are adequately given as a radially hydrostatic equilibrium 
state with zero azimuthal velocities everywhere. Physical 
variables at the inner boundary -Ri n = 2r g , except for the 
velocities, are given by extrapolation of the variables near 
the boundary. However, we impose limited conditions that 
the radial velocities are given by a free-fall velocity and the 
angular velocities are zero. On the rotational axis and the 
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Table 1. Model parameters 



M,/M Q 


Aout 


"out 




Pout (gem 3 ) 


fa <j> 




10 


1.64 


0.0751 


0.0738 


2.1 X 10~ 12 - 7 X 10~- 2 


10~ 4 - 7 x 10 6 29° 


30 


10° 


1.875 


0.0751 


0.0654 


2.1 x 10- 17 - 10~ 6 


10~ 4 - 10 7 29° 


30 



j 1.6E-2 
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Figure 1. Time evolutions of luminosity L (erg s — 1 ) for fa = 
1.6 X 10 _4 ,1.6 X 10 -2 , 16.2, and 7.1 X 10 5 for the stellar-mass 
black hole. 




log T 



8.2 
7.9 
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Figure 3. Same as Fig. 2 but temperature contours. The black 
thick line behind the shocked region shows the funnel wall which 
is characterized by vanishing effective potential. 




Figure 2. Density contours and velocity vectors at t = 1.2 X 
10 4 (r g /c) for fa = 7.1 X 10 5 of the stellar -mass black hole. The 
standing shock is formed at r/r g ~ 8 and z/r g ^3 near to the 
equatorial plane but it is bent obliquely toward upstream. 




Figure 4. Profiles of density p (g cm 3 ), temperature T (K), 
and Mach number of the radial velocity on the equatorial plane 
at t = 1.2x 10 4 (Vg/c) for fa = 7.1 X 10 s of the stellar -mass black 
hole. 



equatorial plane, the meridional tangential velocity w is zero 
and all scalar variables must be symmetric relative to these 
axes. As to the outer boundary region above the outer-disc, 
we use free-floating conditions and allow for outflow of mat- 
ter, whereas any inflow is prohibited here. With these initial 
and boundary conditions, we perform time integration of 
equations (l)-(6) until a quasi-steady solution is obtained. 



4 NUMERICAL RESULTS 
4.1 Stellar-Mass Black Hole 

The input flow with the injection parameters at the outer- 
disc radius R on t arrives at the innermost radius in a time 
smaller than -R ou t/iw(= 40O g /c) and then the system is 
settled down toward a quasi-steady state configuration. The 
luminosity curve is a good measure to check if the quasi- 
steady state is attained. Fig. 1 shows the time evolutions of 
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luminosity L (erg s 1 ) for casese of rh = 1.6 x 10~ 4 ,1.6 x 
10~ 2 , 16.2, and 7.1 x 10 s for the stellar-mass black hole. 
In each case a quasi-steady state is obtained and the cen- 
trifugally supported shock is formed around the black hole. 
However, the luminosities show QPO phenomena with mod- 
ulations of a factor of 2 - 3 and the shock positions on the 
equatorial plane are also found to be variable around several 
Schwarzschild radii. 

The overall features of the flow and the shock are shown 
in Figs. 2 and 3 which show the contours of density p (g 
cm" 3 ) and temperature T (K) with velocity vectors at the 
evolutionary time t — 1.2 x 10 4 (r g /c) for rh — 7.1 x 10 5 , re- 
spectively. Here, a spheroidal shock surface is found, where 
a standing shock exists at r/r g ~ 8 and z/r g ^ 3 near 
to the equatorial plane but it is bent obliquely toward up- 
stream. The densities and the temperatures are enhanced by 
a several factor across the shock surface. Behind the shocked 
region there exists a funnel wall which is characterized by 
vanishing effective potential and is roughly denoted by the 
black thick line in Fig. 3. The extended shocked region be- 
tween the shock surface and the funnel wall consists of the 
hot and dense gas with decelerated velocities. The luminos- 
ity modulations are caused by oscillations of the hot shocked 
region. In the cone-like funnel region between the rotational 
axis and the funnel wall, the temperatures are much higher 
but the densities are very low, and the gas is optically thin. 
The accreting matter in the inner shocked region mostly 
flows into the black hole and partly diverts into the funnel 
region. The outflow gas generated there is originally sub- 
sonic but is accelerated to relativistic velocities due to the 
radiation pressure force. An upper arrow in Fig. 2 shows the 
reference velocity vector of light. 

Fig. 4 denotes the shock profiles of density p (g cm -3 ), 
temperature T (K), and Mach number of the radial veloc- 
ity on the equatorial plane at the same evolutionary time as 
Figs. 2 and 3. The radially infalling gas attains to Mach num- 
ber of ~ 2.7 at the shock and is abruptly decelerated across 
the shock front, and is again supersonically swallowed into 
the black hole. To understand the transitions of flow vari- 
ables at the shock, we refer to the pressure balance equation 
of Rankine-Hugoniot relations: 

Pi + P2V2 2 = -Pi + pwi 2 , (13) 
where P is the total pressure due to radiation pressure £eEi 



fc_ 

and g as pressure R g pT/p, /e is the Eddington factor (iKlevI 
1989), and subscripts "1" and "2" describe quantities before 
and after the shock. 

In the case of rh = 7.1 X 10 , the radiation pressure 
dominates the gas pressure and the gas is optically thick 
everywhere except the funnel region. The shock features in 
Figs. 2-4 are very similar to those in the adiabatic flows 
l|Okuda et al.l 120 04). In Fig. 4, the pre-shock temperature 
of ~ 7.5 x 10 7 K jumps to the post-shock temperature of 
1.5 x 10 s K. Here, / E = 1/3, E = aT 4 , piV > Pi, and 
P2 2> P2V2 2 ■ Then, we have aT2 4 /3 ~ pivi 2 from equation 
(13) and T 2 ~ 1.5 x 10 8 K with Pl = 10" 1 ' 6 g cm -3 and v% = 
0.24. From the numerical data, P2I p\ = vi/vi ~ 5. The gas 
behaves as an adiabatic flow for a perfect gas with 7 = 1.5. 
The flux limiter A expresses the degree of optical thickness of 
the gas and A is 1/3 throughout the whole region except the 
funnel region. On the other hand, for the low accretion rates 
of rh = 1.6 x 10~ 2 and 1.6 x 10~ 4 , the gas is optically thin 




Figure 5. Power energy density spectra of the luminosity L (red 
line) and the shock position R B (blue line) on the equatorial plane 
corresponding to Fig. 1. 




2000 4000 6000 

Time (Rg/c) 

Figure 6. Time evolutions of luminosity L for rh = 1.7 x 10 -4 , 
0.17, 17.0, and 1.7 X 10 7 for the supermassive black hole. 



everywhere and the shock features are considerably differ 
from the adiabatic case, and the shock position near to the 
equatorial plane moves outward compared with those in the 
high accretion rates. 

Fig. 5 denotes the power energy density spectra of lumi- 
nosity L (red line) and shock position i? s (blue line) on the 
equatorial plane corresponding to Fig. 1. From the power 
spectra of L, we recognize the QPO-frequencies of a few to 
10 Hz and the frequencies increase with increasing rh. The 
power spectra of i? s agree qualitatively with that of L. 



4.2 Supermassive Black Hole 

Fig. 6 shows the luminosity curves for casese of rh = 
1.7 x 10" 4 , 0.17, 17.0, and 1.7 x 10 r for the supermassive 
black hole. The luminosities show also the quasi-periodic 
variations with modulations of a factor of 2 - 3. The ab- 
solute luminosities are five orders of magnitude larger than 
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Figure 7. Density contours and velocity vectors at t = 7.4 X 
10 3 (fg/c) for m = 0.17 of the supermassive black hole. The oscil- 
lating shock near to the equatorial plane is formed at r/r g ~ 12 
- 14 and it extends obliquely toward upstream. 




5 10 15 20 25 30 



Figure 8. Same as Fig. 7 but the flux-limiter A which shows the 
degree of optical thickness of the gas. 



those corresponding to the stellar-mass black hole, because 
of the much larger mass of the supermassive black hole. The 
general properties of the flows, the shock profiles, and the 
shock positions for the supermassive black hole are same as 
those for the stellar-mass black hole. The properties of the 
flows depend on the input accretion rates. In the high ac- 
cretion rates, the shock features are same as Figs. 2-4 for 
the stellar- mass black hole and the shock locates at r/r s ~ 
7 - 8 on the equatorial plane, as is thoretically predicted in 
the adiabatic case. However, in the low accretion rates, the 
accretion flows differ considerably from the cases of high m. 

Figs. 7 and 8 show the contours of density and flux- 
limiter A with velocity vectors at t = 7.4 x 10 3 (r g /c) for 
m — 0.17, respectively, where the ambient density p out is as 
low as ~ 10~ 14 g cm -3 . Compared with the flows in Figs. 3 
and 4 for the high accretion rate, the accretion flow is geo- 
metrically thin in the inner region due to the strong radiative 




Figure 9. Profiles of density p (g cm" 1 ), temperature T (K), 
Mach number of the radial velocity, and flux limiter A on the 
eaqatorial plane at t = 7.4 X 10 3 (rg/c) for rh = 0.17 of the 
supermassive black hole. 

cooling effects. The hot and dense shocked region is confined 
in a narrow region (red zone in Fig. 7) above the equatorial 
plane, where the temperatures are as high as 10 8 - 10 10 K. 
Outside of the shocked region, the densities are very low, 
so that the gas is optically thin and the temperatures are 
very high as well. The shock feature near to the equatorial 
plane is found at r/r g ~ 11 - 14 and it extends obliquely 
toward upstream. The optically thick input gas becomes op- 
tically thin (A <C 1/3) in the pre-shock region and is highly 
compressed at the shock, and becomes again optically thick 
(A ~ 0.3) in the post-shock region. 

The profiles of density p (g cm -3 ), temperature T (K), 
Mach number, and flux limiter A on the eqatorial plane are 
shown in Fig. 9. We find here discontinuous structures of 
the flow variables like a shock wave in the range of r/r g = 
11-14 but the flow structures differ considerably from the 
adiabatic shock solution. Since the upstream temperatures 
before the discontinuity are as low as ~ 10 J K due to the 
very low input density, the sound velocity is small and Mach 
number of the infalling gas becomes large. This results in 
large Mach number of ~ 6.4 at the discontinuity. Across the 
discontinuity, the gas is decelerated down to Mach number 
~ 5.2 and is again supersonically falling into the black hole, 
while the infalling gas never become subsonic after the pas- 
sage of the outer sonic point and the flow can not have two 
saddle type sonic points. The discontinuity found in Fig. 9 
is eventually regarded as a pseudo-shock front. However, it 
should be noted that the detailed structures of the horizon- 
tal flows at 5 ^ z /r s ^ 1 above the equatorial plane in Fig. 7 
show the shock features with the two sonic points and that 
the shocked regions are joined to the above discontinuity 
region on the equatorial plane. Then we treat the pseudo- 
shock like a shock wave to estimate the flow variables before 
and after the front. Since the gas pressure dominates the 
radiation pressure in the post-shock region, we have approx- 
imately RgPiTi/ 'fi = pivi 2 . As the result, the low pre-shock 
temperature of ~ 10 5 K jumps to the very high post-shock 
temperature T2 ~ 10 10 K at the shock. The radiation en- 
ergy density never change so drastically as the temperature 
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Figure 10. Power energy density spectra of the luminosity L 
(red line) and the shock position R B (blue line) on the equatorial 
plane corresponding to Fig. 6 

and the density jumps across the shock and is in same order 
of ~ aTi 4 near to the equatorial plane, so that the radia- 
tion pressure dominates the gas pressure in the pre-shock 
region and the gas pressure is dominant in the post-shock 
region. The effective thickness of the shock above the equa- 
torial plane becomes broad instead of the discontinuous one 
in Fig. 4. This agrees with the previous ID result that the 
shock thickness bro adens with decreasing ambient densities 
i|Okudaet al.ll2004h . 

Fig. 10 denotes the power spectra of luminosity L (red 
line) and shock position R s (blue line) on the equatorial 
plane corresponding to Fig. 6. From the power spectra of the 
luminosities, we find the QPO frequencies u qpo = 3 x 10~ 5 
- 3 x 10 -4 Hz and the period -P qpo = 1-10 hours. 



5 SHOCK LOCATION, LUMINOSITY, AND 
MASS-OUTFLOW RATE 

As to the shock waves examined here, we find general prop- 
erties of the shock features independent of the black hole 
masses. This is naturally understood, to some extent, be- 
cause the injection parameters (A ou t, v ou t, flout) have been 
derived from the dimensionless mass, momentum, and en- 
ergy conservation equations for the adiabatic gas which are 
explicitly independent of the black hole masses. There are 
also common properties of normalized luminosities, shock 
positions, and mass-outflow rates versus accretion rate rh 
between the stellar-mass and the supermassive black holes. 
Fig. 11 shows the shock position 7? s versus rh near to the 
equatorial plane for the black holes. For m ^ 10 3 , R B /r g ~ 
7-8, while for rh ^ 1, R B /r g ~ 11 -12. The transition of 
the former shock position to the latter one occurs in the 
range of ro = 1 - 10 2 . The shock positions in the very 
high accretion rates agree with the theoretical adiabatic so- 
lution. In these cases, the gas is optically thick and radiation- 
pressure dominant, and the radiative cooling term balances 
the radiative heating one, that is, the radiative energy source 
A oc (S, — Eo) ~ 0. Accordingly the flow behaves as the adi- 
abatic flow. However, in the low accretion rates where the 



gas is optically thin, S*(— aT 4 ) S> -Eo and the balance of 
cooling and heating is not established anymore. The source 
function A works as the cooling source. At the Rankine- 
Hugoniot relations of a standing radiative shock in the high 
accretion rate, the pressure balance is supported by the sum 
of dominant radiation pressure and ram pressure. When the 
ambient density is taken to be much lower than that in the 
high accretion rate, the upstream pressure just before the 
shock is much smaller because of the lower temperature in 
the pre-shock region. Therefore, to set up a new pressure bal- 
ance condition at the shock, the shock must shift outward 
as far as the same injection parameters are concerened. This 
is the reason why the decreasing input accretion rate leads 
to the increasing shock position on the equatorial plane. If 
the ambient density is taken to be too low, there exists no 
longer shock wave under the same injection parameters, be- 
cause the parameters which are originally derived from the 
adiabatic equations could not match with the optically thin 
flow. 

The QPO-behaviors of the luminosity were attributed 
to the shock oscillations. We furthermore consider that the 
shock oscillations are driven by the centrifug al barrier and 
the r adiative cooling in the shocked region |Molteni et al.l 
1996). When the gas is fully optically thick, the cooling ef- 
fect is negligible. In the transonic accretion problems around 
the black hole, it is well known that the transonic flows have 
generally multi-shock wave solutions, such as the inner and 
outer shock solutions. The shock waves obtained in this pa- 
per correspond to the outer shock wave solution. From the 
instability analyses of accretion flow with a standing shock 
wave, it has been found that the outer shock wave both 
for isothermal and adiabatic flows are dyn amically stable 
(|Nakavamal Il994l ; iNobuta fc Hanawal lilffll ). On the other 
hand, when the flow is in an optically thin state, S* Eo 
and the source function A oc KpT* . As the shock is perturbed 
and propagates outward, it heats the post-shock flow to a 
higher temperature since the relative velocity between the 
shock and the incoming flow becomes higher. As a result, 
the radiative cooling in the post-shock region grows up and 
the outward motion of the shock stops when the flow is suf- 
ficiently cooled down. The shock eventually moves towards 
the black hole due to the reduced pressure behind the shock. 
This time, the relative velocity between the shock and the 
post-shock flow decreases and the post-shock temperature 
drops. The lower pressure in the post-shock region is unable 
to balance the preshock ram pressure and the shock collapse 
continues till the centrifugal barrier is sufficient to hold the 
flow. The shock then bounces outward and this process is 
repeated as the shock oscillations. In order to have an oscil- 
latory behavior, the post-shock region must be able to cool 
in a cooling timescale comparable to the advection timescale 
at the shock. Thus far, the shock oscillation period is con- 
sidered as the advection timescale at the shock position R s 
and the larger R s leads to the larger oscillation timescale 
and the smaller QPO-frequency. From the power spectra in 
Figs. 5 and 10, we find that the QPO-frequencies drift to 
the larger frequencies with increasing rh. 

Fig. 12 shows the normalized luminosity L/Le versus rh 
in the wide range of 10 -4 - 10 7 . The L/Le — rh relations are 
almost identical both for the stellar-mass and the supermas- 
sive black holes. When M is low, the luminosity increases in 
proportion to the accretion rate. However, when M exceeds 
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Figure 11. Shock position i? s on the equatorial plane versus rh Figure 13. Mass-outflow rate Mi obb /Me versus rh for cases of 

for cases of the stellar-mass and the supermassive black holes the stellar-mass and the supermassive black holes. 
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Figure 12. Normalized luminosity L/Le versus rh for cases of 
the stellar-mass and the supermassive black holes. The line shows 
the slim disc solution by Watarai et al. (2000). 

greatly the Eddington critical rate Me, the luminosity is in- 
sensitive to rh and is kept constantly around a maximum 
luminosity of ~ 3Le- We have here approximately 

{log rh - 1.2, for rh sC 10~ 2 
-0.15 (log rh) 2 + 0.74 log m - 0.81, , . 

for 10" 2 < rh < 10 2 , V ; 

0.48, forrh>10 2 . 

IWatarai et all l|200Ch studied the slim-accretion disc 
model shining at the Eddington luminosity and numerically 
derived an L/Le - rh relation. The result is also plotted in 
Fig. 12. Their L/Le - rh relation agrees well with ours in 
the low accretion rates but the luminosities in the high ac- 
cretion rates are rather larger than ours. The result for the 
low accretion rates shows L — 0.063Mc 2 , that is, the con- 
version efficiency r\ of gravitational energy to radiation is 
~ 1/16. The existence of the maximum luminosity ~ 3Le is 
interesting from the ultra luminous X-ray sources point of 
view. The normalized mass-outflow rate mi OS s(= A/i oss /Me) 



versus rh is given in Fig 13. The mass-outflow rate increases 
in proportion to rh and we have 

Afi oss = 0.04 M. (15) 

The mass-outflow rate Mi oss amounts to about a few percent 
of the input accretion rate M. 

6 CONCLUDING REMARKS 

We examined numerically 2D inviscid transonic flows around 
the stellar-mass and the supermassive black holes, while tak- 
ing account of the cooling and heating of gas and radiation 
transport. In these accretion flows, the centrifugally sup- 
ported shock waves are formed around the black holes. The 
shock waves are unstable with time and the resultant lumi- 
nosities show the QPOs with modulations of a factor of 2 

- 3. The shock positions are weakly dependent on the ac- 
cretion rate over a wide range of rh and are in the range 
of 7 - 12 r g on the equatorial plane. In the cases of much 
higher accretion rates (rh 3> 1), the shock position agrees 
with the theoretical adiabatic solution . On the other hand, 
the much lower accretion rate leads to the more outward 
shock position and accordingly the smaller QPO-frequency. 
As results, we have the QPO frequency i/<jpo of a few to 10 
Hz and 3 x 10 -5 - 3 x 10~ 4 Hz, that is, the period P qpo of 
0.1 - several seconds and 1-10 hours for the stellar-mass 
and the supermassive b l ack ho les, respectively. 

IChakrabarti et all l|2004l ) examined numerically the 
shock oscillation model around the black holes with masses 
of IOMq and 1O 8 M0 and obtained the QPO-frequencies of 

- 8 - 12 Hz and 10" 7 - 10" 6 Hz for the stellar -mass and 
the supermassive black holes, respectively. We plot the QPO 
period-mass relation in Fig. 14, together with their results. 
Here we have approximately 

P qpo = 0.016(^) sec. (16) 

This suggests the QPO-periods expected from the shock os- 
cilation model around black holes in various mass range. 
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Figure 14. QPO-period P qpo versus black hole mass M*/M@ 
due to the shock oscillation model. 

The normalized luminosity and the mass-outflow rate 
versus m have common properties independently on the 
black hole masses. When rh is low, the luminosity increases 
in proportion to the accretion rate. However, when M ex- 
ceeds greatly the Eddington critical rate Me, the luminos- 
ity is insensitive to the accretion rate and is kept constantly 
around 3 Le- On the other hand, the mass-outflow rate Mi oss 
increases in proportion to M and it amounts to a few per- 
cent of the input mass-flow rate both for the stellar-mass 
and the supermassive black holes. 

The shock oscillation model works only for the accre- 
tion flows within a limited parameter space with adequate 
input parameters for the inviscid flows. Up to now, only 
the transonic solutions of the accretion flows with total pos- 
itive energy have been discussed and simulated. However, 
recent work of the sub-Keplerian flows with negative total 
energy enlarges the parameter space for which the steady 
shocks are exhib ited in the rotating accretion flows around 
the black holes (|Molteni et al.l [20061 ). Even in the viscous 
accretion discs, the shock oscillation models are still valid 
when the viscosity p a rameter, a, is less than a critical value 
jMolteni et all 1 19961 : lLanzafame et all 119981 1. Of course, if 
the viscosity is too high, the standing shock disappears. 

The centrifugally supported shocks around black holes 
and the shock oscillation models have been applied to 
the QPO p henomena in the b lack hole candidate GRS 
1915+105 JChakrabartl 1 19991 ; IChakrabarti fc Manickaml 
l2000l ; lRao et alJl200Ch in comparison with the observational 
data. The Galactic X-ray transient source GRS 1915+105 
exhibits various types of QPOs: (1) the low-frequency QPO 
(vl ~ 0.001 - 0.01Hz), (2) the intermediate frequency QPO 
(v\ ~ 1 - 10Hz), a nd (3) the high frequency QPO (vn ~ 67 
Hz) (|Morganll 1997m . The QPO-frequencies of a few to 10 Hz 
due to the shock oscillation model in the stellar-mass black 
hole may be representative of the intermediate frequency 
QPO observed in GRS 1915+105. 

The shock oscillation model of the supermassive black 
hole may be applicable to Active Galactic Nuclei, because 
they are powered by accretion onto a central black hole. 
Therefore, it may be plausible that QPOs should be ob- 
served in AGNs as well. The QPO-frequency due to the 



shock oscillation model depends on the mass of the super- 
massive black hole. Recent data of Seyfert galaxies by the 
EUVE and the Chandra show some periodi cities or QPOs 
in the range of a tenth of hours to a month jHalpern et al.l 
l2003tlMoran et alll2005l ) and may suggest the relevance to 
the shock oscillation model. Further observations of the 
black hole candidates will be needed to confirm the shock 
oscillation model for QPOs. 
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